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1. Introduction 



Over the last few years we have performed a study of quenched QCD with overlap quarks. We 
generated 100 configurations of size 18 3 x 64 at /3 = 6 and 100 configurations of size 14 3 x 48 at 
/3 = 5.85, with Wilson gauge action and separated from each other by 10,000 Metropolis steps. The 
corresponding lattice spacings, from the Sommer scale, are a~ l = 2.12GeV and a~ l = 1.61GeV. 
These configurations were transformed to the Landau gauge, and overlap quark propagators were 
calculated for a single point source and all 12 color-spin combinations with p = 1.4, am q = 
0.03,0.04,0.06,0.08,0.1,0.25,0.5,0.75 at J8 = 6, andp = 1.6, am q = 0.03,0.04,0.053,0.08,0.106, 
0.132,0.33,0.66,0.99 at /3 = 5.85. (We recall that the overlap Dirac operator for a quark of mass m 
is given by [1 - (am) / (2p))D + m, where D = (p/a)(l + y*,H (p) / sj H (p) 2 ) and H(p) stands for 
the Hermitian Wilson-Dirac operator with mass —p/a.) A Zolotarev approximation to the inverse 
square root was used for the first 55 configurations at j8 = 6, then a Chebyshev approximation of 
degree 100 ~ 500, after Ritz projection of the 12 (18 3 x 64 lattice) or 40 (14 3 x 48 lattice) lowest 
eigenvectors of H 2 . The convergence criterion was \x~ l l 2 — Y,T n (x) \ < 10~ 8 , \D^D\\f — %\ 2 < 10~ 7 . 
The use of propagators with a point source, as opposed to extended or wall sources, was dictated 
by our desire to calculate matrix elements and renormalization factors, and by the limitation of the 
computational resources at our disposal. Results for light hadron spectroscopy were presented at 
the Lattice 2005 Symposium §\§ and in Ref. [Q]. We refer the reader to [Q] for details of our calcu- 
lations and for references to other studies of lattice QCD with overlap quarks on large lattices. We 
present here results we recently obtained on AS = 2 matrix elements and on baryon wave-functions 
and diquark correlations inside baryons. 



2. AS = 2 matrix elements 

We evaluated matrix elements of the operators 0\ = [j°7^(1 — 7s )<^"] [^7^(1 — ys)d b ],02 = 
- y 5 )d a }{^(\ - y 5 )d b },0 3 = - y 5 )d b ][s b (l - Ys)d a ],0 4 = [s°(l - y 5 )d a ][s b (l + y 5 )d h ] 
and O5 = — y= > )d b ][^ } (1 + ys)d a ] which are needed to study neutral kaon mixing in the stan- 
dard model (SM) and beyond (BSM). An extensive description of our work and results has been 
presented in [^]. Here we we only outline the main features of our analysis. We would like to direct 
the reader to ^ for a detailed bibliography of other relevant investigations. 

We used the quark propagators to calculate 

ail , \ Zx,y(P(x)Oi(0)P(y)) a «v «r/2«vo«r 

Hx,y(P( X )°i(°)P(y)) a «x «r/2«y «T 

NiI,z S (P(x)P(0)){P(0)P(y)) 



?pp(xo,yo) = j. t / p / \^ ^ u B t (2.2) 



for i=2 • • • 5 with A^,- = |, — |, — 2, — |, and fit to a constant in the symmetric time intervals given by: 
12 < x /a < 19 and 45 < y /a < 52 for i=\ ■ ■ -5, at j8 = 6.0; 10 < x /a < 12 and 36 < y /a < 38 
for i=\ and 10 <x /a < 14 and 34 <y /a < 38 for i=2 --5, at j8 =5.85. 

The bare values obtained for the parameters B, must be renormalized to relate them to physical 
observables. We used a non-perturbative renormalization technique based on the RI/MOM methods 
of Our results for the B-parameters are shown in Figure [IJ where the polynomial interpolations 
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Figure 1: Mass-dependence, in terms of the variable M 2 / (4-nF) 2 , of the B-parameters Bi, i = 1, • • • ,5, in the 
RI/MOM scheme at 2 GeV. The solid curves are the results of the fits described in the text, and are plotted 
in the fit region. The fits are used to interpolate the results to the kaon point M 2 / (4-nF) 2 = M\j (\%Fk) 2 , 
shown as a vertical dotted line. The dashed curves are an extension of the fit curves outside the fit range. 
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Figure 2: Same as Figure [j], but for the BSM ratios R l t 



BSM 



to the kaon point are displayed. From these B-parameters, we also reconstruct the matrix elements 
themselves. In Figure ||, we show the polynomial interpolations to the kaon point of the ratios: 







\M 2 {P\Oi{vL)\P\ 






expt 


F 2 (P \O l {lx)\P°)_ 


lot 



(2.3) 



for i = 2, ■ ■ ■ ,5, where M and F are the mass and "decay constant" of the lattice kaon which we 
denote by P° to indicate that the mass of the strange and down quarks that compose it can differ 
from their physical values. The ratios Rf s (n,Mg) measure directly the ratio of BSM to SM 
matrix elements and, as such, can be used in expressions for AM% and e beyond the SM, in which 
the SM contribution is factored out. 

Our main conclusion is that the non-SM, AS = 2 matrix elements are significantly larger than 
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found in the only other dedicated lattice study of these amplitudes [Qj. In tracing the source of 
this difference, we found that we already disagree on the much simpler matrix element of the 
pseudoscalar density between a kaon state and the vacuum, which is the building block for the 
vacuum saturation values of the BSM AS = 2 amplitudes. Through the axial Ward identity, the 
matrix element of the pseudoscalar density is related to the sum of the strange and down quark 
masses, which we find to be roughly 30% smaller than the value obtained in [§], with the same 
tree-level improved Wilson fermion action and gauge configurations as used in [BJ|. Since our 
result for this sum of masses is in agreement with the continuum limit, benchmark result of [R], we 
are convinced that the stronger enhancement of non-SM AS = 2 matrix elements that we observe 
is correct. For details concerning this issue as well as our other results, we refer the reader to [||]. 

3. Baryon wave-functions and diquark correlations 

We study the correlation of quarks inside baryons by evaluating baryon Green functions where 
the quarks, which without loss of generality we take to be u,d and s, are found at different spatial 
locations at the sink: 

G{n,r 2 ,f 3 ,t) = {u(r u t)d(r 2 ,t)s(r 3 ,t) u(0)d(0)s(0)) (3.1) 

The color indices, which like the spin indices are left implicit, are combined in a color singlet. The 
use of gauge fixing allows us to consider separated quarks without the need to introduce gauge 
transport factors. For this investigation we also converted the gauge background field from the 
Landau gauge to the Coulomb gauge, and we will report below results in both gauges, which 
however are rather similar. In order to project over states of zero momentum, we sum over a 
translational degree of freedom, evaluating a reduced wave function 

G(r/,t)='£{u(r 3 + r,t)d(r 3 +r',t)s{r 3 ,t) m(0)</(0)j(0)) (3.2) 



(The sums in involve a very large number of terms and the use of a fast Fourier transform and 
the convolution theorem is crucial to carry them out in manageable computer time). 

The spin indices are combined in an appropriate spin wave-function at the source and sink. 
Of particular interest is comparing the spin configurations where, in the spin l/2 + octet, the u 
and d quarks are combined, A-like , in a spin and isospin singlet state, the so-called good diquark 
combination, versus those where u and d are, Z-like, in a spin and isospin triplet state. In Figure || 
we show the ratio of the mean separations between u,d and u,s quarks in the two spin states. The 
results give support to the notion that quarks inside a baryon tend to correlate in a "good" diquark 
state. 

It is interesting to visualize the correlation among the quarks inside the baryon, i.e. the function 
G(r,r',t) which we can loosely consider as the wave-function of the three quarks. The problem 
is, of course, that, at fixed t, G is a function of the two vectors r,r', i.e. of six variables. The 
representation of G can be simplified somewhat by assuming rotational invariance, i.e. lack of 
major spin orbit correlation. We have verified that this assumption is satisfied within the statistical 
accuracy of our calculations. (Indeed the actual magnitude of spin-orbit correlations could be 
evaluated by our technique, given sufficient statistical precision.) This leaves G a function only of 
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Figure 3: Ratio of mean u — d separation to mean u — s separation as function of quark mass for the two 
u,d diquark configurations. Data with j3 = 6 and t = 10a. 



Figure 4: The geometry used for visualizing the wave-function of the quarks inside a baryon. 



the shape of the triangle subtended by the two vectors f, f . This is still a function of three variables, 
but, to provide a meaningful visualization, we can fix one of these and represent the wave-function 
as a function of the other two. Thus, for a generic triangle formed by the locations of the s, u and d 
quarks, we introduce coordinates x,y,z as illustrated in Figure |], and then represent G as a function 
in the x,y plane at fixed z. Our results, with /3 = 6, am q = 0.03 and t = 10a, are illustrated in 
Figures || and []. They show again that the u and d quarks tend to correlate in the good diquark 
configuration. 

A detailed, expanded version of the results presented in this section is in preparation and will 
form the subject of a forthcoming paper. An earlier study of quark wave-functions with some 
similarity to ours can be found in Ref. [Qj. We are not aware of other investigations studying the 
same type of correlation functions. 
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Figure 5: The wave-function of the u and d quarks, in the Coulomb gauge, in a spin and isospin singlet 
state (good diquark), given by the narrower bell-shaped surface (red with color), against the wave-function 
of the two quarks in a spin and isospin triplet state, given by the broader curve (green with color), for a 
separation z = 1 -8a between the s quark and the mid-point of the u, d pair. 



4. Conclusions 

Our results corroborate the fact that the overlap discretization, at least insofar as valence quarks 
are concerned, can be used in large scale simulations, and, because of its very good symmetry 
properties, represents a choice method for QCD numerical calculations. 

They also provide some novel matrix element values and evidence for strong diquark correla- 
tion in a flavor 3, spin-singlet state. 
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Figure 6: Same as in Figure 5, but for z = 3.6a. 
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